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INTRODUCTION 


When a laminate containing a near-surface delamination is subjected to compression 
loads, local buckling of the delaminated region can occur. This buckling causes load re- 
distribution and secondary loads, which in turn cause interlaminar stresses. Delamination 
growth resulting from these interlaminar stresses is referred to herein as instability-related 
delamination growth (IRDG). Fig. 1 shows two configurations which exhibit IRDG: uniax- 
ially loaded laminates with a through-width or embedded delamination. Typical analytical 
studies of these configurations can be found in refs. 1-8. The case of a circular embedded 
delamination with axisymmetric loads has also been examined (ref. 9). The mechanics of 
IRDG have been described for the through-width delamination in refs. 2 and 4. Herein, 
the term “mechanics of IRDG” refers to the process by which load is redistributed and 
secondary loads are created. A description of the mechanics of IRDG includes much more 
than assemb lying the governing equations, solving them, and obtaining some values for 
the strain-energy release rates. Unless an intuitive feel for IRDG is acquired, even exact 
answers for idealized configurations are of limited use. Earlier studies of the embedded 
delamination have presented various analyses, but the mechanics of IRDG have not been 
addressed. The purpose of this paper is to describe the mechanics of IRDG for the em- 
bedded delamination. 

Because the stresses are singular along the delamination front, point values of the 
stresses are of limited use. Strain-energy release rates are finite parameters which 
quantify the intensity of the stresses at the delamination front. In this paper it will 
be assumed that the magnitude of the strain-energy release rates govern delamination 
growth. Consequently, the mechanics of IRDG will be discussed in terms of its influence 
on strain-energy release rates. 

The stress analysis used in this study was a geometrically nonlinear three-dimensional 
finite element program (NONLIN3D, ref. 8). The following sections will begin with a 
brief description of NONLIN3D, the finite element models, and the material properties. 
Then typical strain-energy release rate results will be shown for four configurations 
which exhibit IRDG. The four configurations are the through-width delamination, the 
embedded delamination with axisymmetric loads, and circular and elliptical embedded 
delaminations with uniaxial loads. These strain-energy release rates will be used to 
illustrate the distinct differences in behavior exhibited by the configurations. Then the 
mechanics of IRDG will be discussed. First the mechanics of IRDG for the through-width 
delamination will be reviewed. Then its close relationship to the axisymmetrically loaded 
embedded delamination is explained. A laminate with a through-width delamination can 
be transformed into a laminate with an embedded delamination by applying appropriate 
tractions. The mechanics of IRDG for the embedded delamination will be explained by 
considering the type and sign of traction required. The differences in the strain-energy 
release rates for the four configurations will be discussed in terms of the mechanics of 
IRDG. 


NOMENCLATURE 


a, b 

Semi-axes of elliptical delamination in x and y 
directions, respectively 

c a 

Constitutive coefficients 

#lli #22> #33 

Young’s moduli for orthotropic material 

Fy,Fz 

Forces in y- and z-directions, respectively 

G i, G n 

Mode I and mode II total strain-energy release rate 

G Vh G 23j G 13 

Shear moduli for orthotropic material 

h 

Thickness of sublaminate 

H 

Thickness of base laminate 

My 

Moment about y axis 

U, V, w 

Displacements in x-, y-, and z-directions 

W 

Length and width of the square finite element model 

x, y, z 

Cartesian coordinates 

a 

Angle between local and global coordinate systems 

S 0 

specified axial strain 

* / 12» ^23> ^13 

Poisson’s ratios for orthotropic material 
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ANALYSIS 


This section will give a brief description of the finite element program used to obtain the 
strain-energy release rates presented herein. The finite element program, which is named 
NONLIN3D was developed at NASA Langley Research Center. Further details may be 
found in ref. 8. Also, a typical finite element model and the material properties will be 
described. 

The Finite Element Program NONLIN3D 

NONLIN3D is a three-dimensional, geometrically nonlinear finite element program. The 
primary type of geometric nonlinearity considered is that due to significant rotations. This 
is accounted for by using the Lagrangian nonlinear strain-displacement relations (ref. 10). 
The governing equations were derived by minimization of the total potential energy. A 
Newton-Raphson iterative procedure was used to solve the nonlinear governing equations. 

Figure 2 shows two typical finite element models. The same 2D mesh was used to 
generate both of the models shown. Hence, the close-up in fig. 2a is also valid for fig. 2b. 
The elements are 20-node isoparametric hexahedra. The sublaminate thickness h was 
one-tenth of the base laminate thickness H. Because of symmetry only one-fourth of the 
specimen was modeled and the constraints u = 0 on x = 0 and v = 0 on y = 0 were 
imposed. There was also a constraint w = 0 on 2 = 0. This constraint was imposed 
to remove global bending from the analysis. In reality, there might be global bending 
(particularly if the buckled region is thick), but the amount of global bending would 
depend on the size of the region modeled and the boundary conditions at the external 
boundaries. Imposing w = 0 on z — 0 simply removed overall specimen size and external 
boundary conditions as parameters to be considered in this study. The constraint on w 
represents a laminate which is well constrained globally. Of course, one could also view 
the constraint u; = 0onz = 0asan indication of symmetry about the z = 0 plane. This 
implies the presence of two delaminations. 

Along the boundary x = W, all u displacements are specified to equal xe 0 , where e 0 
is the specified compressive axial strain. To initiate transverse deflections, a transverse 
load was applied at the center of the delaminated region. After a converged solution 
was obtained, the load was removed, and solutions were obtained with only compression 
loading. 

The model for the axisymmetric load case differed only in that equal displacements were 
prescribed along both x=W and y=W to obtain biaxial compression. This biaxial loading 
results in axisymmetric loading of the delaminated region. 

The model for the through-width delamination is shown in fig. 2b. Since this case is 
basically 2-D, a 2D analysis was performed. In order to perform a two-dimensional analysis 
using NONLIN3D, v was prescribed to be zero on the planes y=0 and y=W. This restricted 
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the model behavior to 2-D plane strain. The through-width delamination case is trickier 
to analyze than the other cases because the tangential stiffness of the buckled sublaminate 
is essentially zero. This results in a singular stiffness matrix. Refs. 2 and 11 describe the 
procedure used to circumvent this problem. 

To reduce the computational effort, substructuring was used. Most of the postbuckled 
region was in one substructure and the remainder of the laminate was in the other 
substructure. Fig. 3 shows a typical division of the laminate. The boundary between 
the two substructures was a distance “£” behind the delamination front. The length £ 
was approximately equal to the sublaminate thickness h . Since significant rotations occur 
only in the buckled region, a nonlinear analysis was required only in one substructure. 
A reduced stiffness matrix and load vector were obtained for the linear substructure and 
these were used to augment the equations for the nonlinear substructure. 

Strain-energy release rates were the primary output from the analysis. The strain- 
energy release rates were calculated using a three-dimensional version of the virtual crack 
closure technique in ref. 13. The strain-energy release rates were calculated relative to 
a local coordinate system. The local coordinate system for the circular and elliptical 
delaminations is shown in fig. 4. One axis is normal to the delamination front and one axis 
is tangent. For all the cases considered, the local z-axis was parallel to the global z-axis. 
The angle between the local and global y-axes is a. 

Material Properties 

The intent of this paper is to consider only the effect of geometric parameters on the 
mechanics of IRDG. Consequently, material properties were chosen to minimize the effect 
of material properties. For quasi-isotropic laminates the in-plane stiffness is independent of 
direction. But even for quasi-isotropic laminates the flexural stiffness varies with direction. 
Hence, even if the postbuckled region consisted of a quasi-isotropic group of plies, one would 
expect variations in the strain-energy release rate along the delamination front which are 
caused solely by the variation in flexural stiffness. Also, the properties of the interface 
plies (i.e., those plies on either side of the delamination) would be expected to at least 
affect the relative magnitudes of Gj, Gjj, and Gjjj. 

The simplified material properties chosen for this study are those for a "homogeneous 
quasi-isotropic laminate”. These properties C t y are obtained as follows: 

= i 

t=l 


where [C^jY are the constitutive properties for the £ 4 ^ ply in the 8-ply quasi-isotropic 
laminate (±45/0/90) s . With these properties, there are no stacking sequence effects and 
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no variation of material properties with orientation. These properties are not isotropic, 
since there are still reduced transverse shear and extension moduli. 

The lamina properties were selected to be typical of graphite/epoxy (ref. 12 ). The 
moduli and Poisson’s ratios are: 


E\\ = 134 GPa 
E 22 — £33 = 10.2 GPa 
G n = G 13 = 5.52 GPa 
C ?23 = 3.43 GPa 
u 12 = ^13 = -3 
^23 = -49 

COMPARISON OF BEHAVIORS FOR THE THROUGH- WIDTH AND 
THE EMBEDDED DELAMINATION 


Four cases were analyzed: the through-width delamination, the axisymmetrically loaded 
embedded circular delamination, the uniaxially loaded embedded circular delamination, 
and the uniaxially loaded embedded elliptical delamination. In all cases the thicknesses 
H and h were 4 and .4 mm, respectively. The through-width delamination was 30mm 
long. The circular delaminations were 30mm in diameter. The elliptical delamination was 
elongated along the y-axis and had dimensions of 30x60mm. The same strain range was 
used for all the cases. 

Figures 5 and 6 show the effect of strain level on Gj and Gjj for the through-width 
delamination and the embedded delamination. The mode III component is not shown 
because it was essentially zero for all the cases considered. In general, one would not 
necessarily expect Gjjj to be zero. Both uniaxial and axisymmetric loads were considered 
for the embedded delamination. For the uniaxially loaded embedded delamination both 
a circular and an elliptical delamination were analyzed. Figs. 5 c, 5 d, 6 c, and 6 d show 
the variation with strain level at two points along the delamination front of the uniaxially 
loaded embedded delamination: at a = 0 ° and a — 90°. See fig. 4 for the definition of 
alpha. 

The variation of Gj with strain level is dramatically different for the four cases. For the 
through-width delamination Gj increases very rapidly after buckling occurs. (fig. 5 a) After 
reaching a peak, Gj decreases. For the axisymmetric case, Gj increases monotonically. 
Also, the magnitude of Gj in fig. 5b is much larger than in fig. 5 a. 

Fig. 5c shows the variation of Gj at a = 90° for a uniaxially loaded embedded 
delamination. Only the results for the elliptical delamination are shown here, since Gj 
wa s zero at a = 90° for all strain levels for the circular delamination. The shape of the 
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curve for the elliptical delamination is similar to that for the through-width delamination 
(fig. 5a), but the magnitude is less. At a = 0°, Gj is shown for both the circular and 
elliptical delaminations (fig. 5d). Note that the magnitude of Gj is much larger than in 
figs. 5a-5c. Also, Gj increases rapidly and monotonically with applied strain. 

Fig. 6 shows the Gjj variations with strain level. In all cases G// increased 
monotonically with strain. G// is of the same order of magnitude for all the cases except 
for the uniaxially loaded circular delamination at a = 90° (fig. 6c). This contrasts with 
the very wide range of magnitudes in fig. 5 for Gj. 

Even though the strain-energy release rates in all the cases illustrated in figs. 5 and 6 
are a result of local buckling, the variety of behavior suggests that there must be variations 
in the mechanism by which local buckling causes strain- energy release rates. The next 
section will attempt to explain these mechanisms. 

MECHANICS OF IRDG FOR THE EMBEDDED DELAMINATION 

In highly simplified anthropomorphic terms, a strip of the buckled region which is 
parallel to the load direction (strip A in fig. 7) wants to buckle outward. A strip of the 
buckled region which is perpendicular to the load direction (strip B in fig. 7) has no desire 
to deform outward; it is pushed out by strip A. Strip A is analogous to the through-width 
case. The constraint provided by strip B reduces G/ for strip A. Conversely, strip A 
causes high Gj at the ends of strip B when strip A pushes strip B outward. Of course, the 
buckled region is not comprised of strips, but this simplified interpretation helps explain 
the behavior observed. The following paragraphs present a more rigorous and detailed 
discussion of the mechanics of IRDG. 

The through-width delamination will be discussed first. After describing the mechanics 
for the through-width delamination, its close relationship with the embedded delamination 
with axisymrnetric loads will be discussed. Next, tractions will be applied to the through- 
width delamination configuration which transform it into a uniaxially loaded embedded 
delamination. The required tractions should give some feel for why the behaviors differ 
for the embedded delamination and the through-width delamination under uniaxial loads. 

The discussion of the through-width delamination can be expedited by first transforming 
this geometrically nonlinear problem into a linear one with nonlinearly related loads 
(ref. 4). Fig. 8a shows a schematic of a laminate with a postbuckled through-width 
delamination. In fig. 8b the buckled region is replaced by the loads Pjj and M, the axial 
load and moment, respectively, in the buckled region where it is cut. The total applied 
load Pj> is equal to Pq + Pq • The load system in fig. 8c, (which is the same as fig. 8b) can 
be divided into the two load systems shown in figs. 8d and 8e. Because Pq and Pg are 
calculated using rule of mixtures, the load system in fig. 8e causes a uniform axial strain 
state and no interlaminar stresses. Accordingly, for strain-energy release rate calculation, 
the configuration in fig. 8d is equivalent to the original configuration (fig. 8a). The moment 
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M opens the delamination, contributing to Gj. It also contributes to some Gjj (ref. 4). 
The load ( P c - Pd) contributes to Gjj. In addition, because of the offset of the line 
of action of Pq — Pd relative to the delamination, this force creates a moment which 
tends to close the delamination and reduce the Gj component caused by M. The result of 
the competing mechanisms are strain-energy release rate variations like that in fig. 5a and 
6a. The mode I strain-energy release rate first increases very rapidly with increasing strain 
and then decreased to zero. The mode II strain-energy release rate increases monotonically 
with applied load, since both M and ( P c - P D ) contribute to G n . 

The axisymmetrically loaded circular delamination is very similar to the through-width 
delamination. In fact, the schematics in fig. 8 are applicable if the forces are replaced by 
forces per unit length. The load in a column, P p is essentially constant after the applied 
strain is increased beyond the buckling strain, but the load in an axisymetrically loaded 
plate continues to increase significantly after buckling. The load Pq increases linearly 
with the applied load. Hence, {Pq-Pd) is large for the through-width delamination but is 
relatively small for the axisymmetric case. As a result there is little attentuation of the 
effects of M by ( P c - P D ) for the axisymmetric case. Consequently, Gj is much larger for 
the axisymmetric case (Compare figs. 5a and 5b). 

Now the uniaxially loaded embedded delamination will be considered. Figs. 9-12 
illustrate the transformation of a through-width delamination (fig. 9) into an embedded 
delamination (fig. 12). The letters A through G are added to aid the discussion. They are 
not related to the letters in fig. 8. To expedite the discussion, the embedded delamination 
will be assumed to be rectangular. Tractions are require to close the buckled part of region 
AEFB. These will only be non-zero near the new delamination front and are, in fact, the 
interlaminar stresses. 

Fig. 10 shows a slice removed from the laminate. The figure shows the forces required 
to close the buckled part of the boundary BF. These forces are generated when the region 
AEFB is closed. These forces indicate some of the interaction of regions AEFB with BFGC. 
There are in-plane forces F y , transverse forces F z , and a moment M y . Fig. 11 shows the 
same slice after the forces have closed the buckled part of BF. 

The moment M y , would operate in the direction indicated in fig. 10 based on the 
curvature in fig. 11. Likewise, the transverse force, F z , would be expected to act downward 
to help close the delamination front. The force, F y , is a result of two things : transverse 
deflection and Poisson’s ratio. When transverse deflection occurs, the length of a line from 
A to D (fig. 12) increases, hence the buckled region must be stretched in the y-direction. 
When the laminate is compressed in the x- direction, it expands in the y-direction due to 
Poisson’s effect. If the base laminate has a larger Poisson’s ratio than the sublaminate, 
then a force F y is required to enforce compatibility when the buckled part of BF is closed. 
The sign of F y due to Poisson’s ratio would depend on the relative magnitudes of the 
Poisson’s ratios. Differences in Poisson’s ratio were not considered in this study. 
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The magnitude of F y is should be related to the in-plane stiffness in the y-direction. The 
magnitude of My should be related to the flexural stiffness in the y-direction. The effect 
of material properties on F z is not as straight-forward, so no prediction will be offered. 

The dimensions of the embedded delamination should affect the forces and moment. 
For the same transverse deflection, the curvature in the y-direction is less for a larger 
dimension b , so M y should decrease with increasing b. Also, the strain in the y-direction 
due to the increased length A-D (fig. 12) would be less for larger b. Hence, for the same 
transverse deflection F y should decrease with increasing b. 

The moment M y should contribute primarily to Gj, but it also contributes to G/j. 
The force F y should contribute to G// and reduce G/. The reduction is due to the offset 
between the delamination plane and the middle of the sublaminate. This offset causes a 
moment relative to the delamination plane which is opposite to M y . Based on the large 
Gj in fig. 5d, the opening effects of M y must dominate the closing effects of F y . 

The original configuration with a through-width delamination (fig. 9) had some distri- 
bution of Gj and G// along x= a. The application of the forces F y , F z , and the moment 
M y changes the load flow significantly. Closing the ends of the through-width delamination 
(area AEFB) contributes a compressive component of interlaminar stresses along the front 
x=a, thus reducing G/. In fact, because of this reduction, G/ was zero for all strain levels 
at a — 90 0 for the circular delamination under uniaxial loads. For the elliptical delamina- 
tion, the magnitude of G/ at a = 90° was less than for the through-width delamination 
for the same reason (fig. 5a and 5c). 

Based on the preceeding discussion, one would expect the behavior of a through-width 
delamination and the uniaxially loaded embedded delamination to have some similarity, 
but probably more differences. Fig. 5 illustrates this very well. For an embedded 
delamination highly elongated perpendicular to the load direction, one would expect 
the behavior to be like that for the through-width delamination. No highly elongated 
delaminations were examined in this study, but even the 1:2 aspect ratio ellipse has a Gj 
variation at a = 90° (fig. 5c) which is very similar to that for a through-width delamination 
(fig. 5a) . But this similarity has little importance for this case, since the Gj was very 
much larger at a = 0° (fig. 5d). 

SUMMARY 

The variation of strain-energy release rate with applied strain was shown to be 
dramatically different for the through-width delamination and the embedded delamination. 
Both axisymmetric and uniaxial loads were considered for the embedded delamination. 
These divergent behaviours were explained in terms of the mechanisms which create the 
strain-energy release rates. 

Earlier studies described the mechanics of instability-related delamination growth 
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(IRDG) for the through-width delamination. For completeness these earlier results were 
discussed herein. This paper also describes the mechanics of IRDG for the embedded 
delamination. The mechanics for the axisymmetric case were found to be very similar to 
that for the through-width delamination. Large differences in Gj observed for the through- 
width and axisymmetric cases could be explained by the different mix of secondary forces 
that generate strain-energy release rate. The uniaxially loaded embedded delamination 
has a much more complicated response than either the through-width or axisymmetric 
cases. By starting with a through-width delamination and applying tractions, one can 
transform it into an embedded delamination. Examination of these tractions explains 
why the uniaxially loaded embedded delamination behaves very differently than either the 
through-width or the axisymmetric cases. 
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Fig. 1 Two configurations which exhibit instability-related 
delamination growth. 







(a) Embedded delamination 


Fig. 2 Typical finite element models. 



delamination front 


Fig. 3 Typical model after division into substructures. 



(a) Through-width delamination 


(b) Embedded delamination with 
axisymmetric load 



(c) Embedded delamination with 
uniaxial load 

Fig. 5 Effect of strain level on strain-energy 



(d) Embedded delamination with 
uniaxial load 

release rate 




(a) Through-width delamination 



(c) Embedded delamination with 
uniaxial load 



(b) Embedded delamination with 
axisymmetric load 



(d) Embedded delamination with 
uniaxial load 


Fig. 6 Effect of strain level on mode II strain-energy release rate 
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Fig. 7 Strips of sublaminate aligned perpendicular and parallel 
to the load direction. 
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(a) 


(b) 



(c) 


(d) 


(e) 


Fig. 8 Nonlinear configuration (a) transformed into linear configuration 
(d) with two nonlinearly related loads (Pc~ Pp) and M. 
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Fig. 9 Schematic for transformation of through-width delamination 
into an embedded delamination. 
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Fig. 10 Slice of laminate showing tractions required to perform 
transformation by closing buckled part of boundary BF 
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F 


Fig. 11 Slice of laminate after tractions have closed boundary 
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Fig. 12 Entire laminate after transformation to embedded delamination 
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